^Explanation of the similarity of the photoemission spectra of SrVOs and CaVC>3 
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We present parameter-free LDA+DMFT results for the many-particle density of states of cubic 
SrVC>3 and orthorhombic CaVC>3 . Both systems are found to be strongly correlated metals, but not 
on the verge of a metal-insulator transition. In spite of the considerably smaller V-O-V bonding 
angle in CaVOa the photoemission spectra of the two systems are very similar, their quasiparticle 
parts being almost identical. This is in contrast to earlier theoretical and experimental conclusions, 
but in agreement with recent bulk-sensitive photoemission experiments. 
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Transition metal oxides are an ideal laboratory for the 
study of electronic correlations in solids. Among these 
materials, cubic perovskites have the simplest crystal 
structure and thus may be viewed as a starting point 
for understanding the electronic properties of more com- 
plex systems. Typically, the 3d states in those materials 
form comparatively narrow bands of width 2 — 3 eV 
which lead to strong Coulomb correlations between the 
electrons. Particularly simple are transition metal ox- 
ides with a 3d 1 configuration since they do not show a 
complicated multiplet structure. 

Intensive experimental investigations of spectral and 
transport properties of strongly correlated 3d 1 transi- 
tion metal oxides started with the paper by Fujimori et 
al. . These authors observed a pronounced lower Hub- 
bard band in the photoemission spectra (PES) which can- 
not be explained by conventional band structure theory. 
A number of papers subsequently addressed the 
spectral, transport and thermodynamic properties of the 
3d 1 series Sri-zCa^VC^ for various values of x, yielding 
contradictory results. While the thermodynamic prop- 
erties (Sommerfeld coefficient, resistivity, and paramag- 
netic susceptibility) were found to be essentially inde- 
pendent of x from x = (SrV0 3 ) to x= 1 (CaV0 3 ) 0, 
PES and Bremsstrahlungs isochromat spectra (BIS) 
showed drastic differences. In fact, the spectroscopic 
data seemed to imply that this substitution series de- 
velops from a strongly correlated metal (SrVOa) into, 
practically, an insulator (CaVOa), i.e., that for x — > 1 
Sri_2;Ca x V03 is on the verge of a Mott-Hubbard transi- 
tion. 

An experimental answer to this puzzle was very 
recently provided by bulk-sensitive PES obtained by 
Maiti et al. 0, and especially by high-resolution bulk- 
sensitive PES measured by Sekiyama et al. 0. In the 
latter work it was shown that (i) the technique of prepar- 



ing the sample surface (which should preferably be done 
by fracturing) is very important, and that (ii) the energy 
of the X-ray beam should be large enough to increase the 
photoelectron escape depth to achieve bulk-sensitivity; 
also the beam should provide a high instrumental reso- 
lution (about 100 meV in 0). With these experimental 
improvements the PES of SrV03 and CaV03 were found 
to be almost identical 0,0, implying consistency of spec- 
troscopic and thermodynamic results at last. This is also 
in accord with earlier Is x-ray absorption spectra (XAS) 
by Inoue et al. which differ only slightly above the 
Fermi energy, in contrast to the BIS data 0. 

The main effect of substituting Sr ions by the isovalent, 
but smaller, Ca ions is to decrease the V-O-V angle from 
(9 = 180° in SrVOg to (9«162° in the orthorhombi- 
cally distorted Pbnm structure of CaV03 . This bond 
bending results in a decrease of the one-particle band- 
width W and thus in an increase of the ratio U/W as 
one moves from SrVC>3 to CaV03. Within a one-band 
Hubbard model, neglecting the orbital structure of the V 
3d, Rozenberg et al. 01 were able to fit the more surface- 
sensitive spectroscopic data for Sri-zCa^VC^ and, with 
other parameters, the bulk-sensitive PES 0. 

In this Letter we present results of a parameter- 
free, comparative study of SrV03 and CaV03 us- 
ing the recently developed LDA+DMFT computational 
scheme B 13 11 E 0, which merges the 
local density approximation (LDA) with a modern 
many-body technique, the dynamical mean-field theory 
(DMFT) [rj| ■ In this realistic approach the orbital de- 
grees of freedom of the three partially filled t2 9 bands 
are explicitly taken into account 0, resulting in a spec- 
tral weight distribution distinctly different from that of 
a one-band model. 

We first discuss the electronic and crystal structure of 
SrV03 and CaVC-3. The valence states of these two sys- 
tems consist of completely occupied oxygen 2p bands and 
partially occupied V-3d bands. There is one electron oc- 
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FIG. 1: Partial DOS of the V-3d t 2g band (solid line) and e 9 
band (dashed line) for SrVOs and CaVOs calculated by the 
LDA-LMTO method. The Fermi level corresponds to zero 
energy. 

cupying the d-states per V ion (d 1 configuration). Fig. Q] 
shows the LDA density of states (DOS) for both com- 
pounds which we obtained via the tblmto47 code of 
Andersen and coworkers 0, 0] , using the crystal struc- 
ture data of Refs. 0,0. 

In both compounds the V ions are in a octahedral oxy- 
gen coordination. The octahedral crystal held splits the 
V-3d states into three degenerate t 2g and two degenerate 
e g states. While a hybridization between these t 2g and c g 
states is forbidden for the cubic SrV0 3 crystal, the dis- 
torted orthorhombic structure of CaV03 allows them to 
mix. In Fig. ^ one can see some contributions of the t 2g 
and e g subbands in the energy range from -8 eV to -2 eV 
due to the hybridization with 0-2p states. They amount 
to 12% and 15% of all t 2g states in SrV0 3 and CaV0 3 , 
respectively. The t 2g states are dominant in the vicinity 
of the Fermi energy, the center of gravity of the e g band 
lying above the upper band edge of the t 2g band. There- 
fore the low energy physics (< 1.5 eV) is governed by the 
t 2g bands which can be considered sufficiently separated 
from the e g states in both compounds. 

We concentrate on this energy range in Fig. [21 where 
we compare the LDA t 2g DOS of SrV0 3 and CaV0 3 . 
Most importantly, the one-electron t 2g bandwidth of 
CaV0 3 , defined as the energy interval where the DOS 
in Fig.^is non-zero, is found to be only 4% smaller than 
that of SrV0 3 (W Ca vo 3 = 2.5 eV, W SrV o 3 = 2.6 eV). 
One would have expected that the strong lattice distor- 
tion with a decrease of the V-O-V bond angle from 180° 
to 162° affects the t 2g bandwidth much more strongly. 
Such a larger effect indeed occurs in the e g bands whose 
bandwidth is reduced by 10%. To physically understand 
the smallness of the narrowing of the t 2g bands we have 
calculated the effective t 2g -t 2g and e g -e g hopping parame- 
ters. The predominant contribution to the e g -e g hopping 
is through a d—p—d hybridization, which is considerably 
decreasing with the lattice distortion. This is also the 
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FIG. 2: Comparison of the LDA DOS of the V-3d t 2g band 
for SrVOs (solid line) and CaVOs (dashed line). In CaVOs 
the degeneracy of the three t2 a bands is lifted (see inset). 
The bandwidth WcaV0 3 = 2.5eV is only 4% smaller than 
W S rvo 3 = 2.6 eV. 

case for the t 2g orbitals. However, for the t 2g orbitals 
the direct d—d hybridization is also important. This hy- 
bridization increases with the distortion since the t 2g or- 
bital lobes point more directly towards each other in the 
distorted structure. Altogether, the competition between 
i) decreasing d—p—d hybridization and ii) increasing d—d 
hybridization results in a very small change of the t 2g 
bandwidth. This explains why previous suggestions of 
strongly different bandwidths are untenable. 

As is well-known, LDA does not treat the effects of 
strong local Coulomb correlations adequately. To over- 
come this drawback, we use LDA+DMFT as a non- 
perturbative approach to study strongly correlated sys- 
tems ^3 EH- It combines the strength of the LDA in 
describing weakly correlated electrons in the s- and p- 
orbitals, with the DMFT treatment of the dynamics due 
to local Coulomb interactions. In the present paper we 
will discuss the relevant parts of the LDA+DMFT ap- 
proach only briefly, refering the reader to Ref. [Hi ] for 
details. 

Given a specific material, it is possible to extract 
from the LDA band structure a one-particle Hamilto- 
nian -H^da wnere the averaged Coulomb interaction is 
subtracted to avoid double counting ^l| ■ Supplementing 
-Hlda with the local Coulomb interactions between the 
electrons one arrives at a material-specific Hamiltonian 
which includes correlations: 

H = HluA + ^EE (i) 

m i 

~T" ^ ^ ^ ^ ^ {U fiaa' J) flima^iva' a' ■ 

In the present case, the index i enumerates the V sites, 
m denotes the individual t 2g orbitals, and a is the spin. 
Because of the nearly cubic symmetry of CaV0 3 , one 
can simplify the calculation and use only the one-particle 
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FIG. 3: LDA+DMFT(QMC) spectrum of SrVOa (solid line) 
and CaVC>3 (dashed line) calculated at T=300 K (inset: effect 
of temperature in the case of CaVOs). 

LDA DOS N°(e) of Fig-Hinstead of the full Hamiltonian 
-^lda El- ^ n Eq- CJ, the local intra-orbital Coulomb 
repulsion [/, the inter-orbital repulsion U', and the ex- 
change interaction J are explicitly taken into account. 
We calculated these interaction strengths by means of 
the constrained LDA method [2(| for SrV03, allowing the 
c g states to participate in screening j^J. The resulting 
value of the averaged Coulomb interaction is U = 3.55 eV 
(U = U' for t 2g orbitals [3, H3) and J = 1.0 eV. The 
intra-orbital Coulomb repulsion U is then fixed by rota- 
tional invariance to U = U' + 2 J = 5.55 eV. We did not 
calculate U for CaVC-3 because the standard procedure 
to calculate the Coulomb interaction parameter between 
two t2 g electrons which arc screened by e ff states is not 
applicable for the distorted crystal structure where the 
c g and t2 g orbitals are not separated by symmetry. On 
the other hand, it is well-known that the change of the 
local Coulomb interaction is typically much smaller than 
the change in the DOS, which was found to depend only 
very weakly on the bond angle. That means that U for 
CaV03 should be nearly the same as for SrV03. There- 
fore we used U = 3.55 eV and J = 1.0 cV for both SrV0 3 
and CaV03. This is also in agreement with previous cal- 
culations of vanadium systems |2l| and experiments [2(| . 

The Hamiltonian is solved within the DMFT us- 
ing standard quantum Monte-Carlo (QMC) techniques to 
solve the self-consistency equations |2jJ • From the imagi- 
nary time QMC Green function we calculate the physical 
(real frequency) spectral function with the maximum en- 
tropy method. The resulting LDA+DMFT(QMC) spec- 
tra for SrVOs and CaV03 in Fig. show genuine corre- 
lation effects, i.e., the formation of lower Hubbard bands 
at about —1.5 eV and upper Hubbard bands at about 
2.5 eV, with well pronounced quasiparticle peaks at the 
Fermi energy. Therefore both SrV03 and CaV03 are 
strongly correlated metals. The 4% difference in the LDA 
bandwidth between SrV03 and CaV03 is only reflected 
in some additional transfer of spectral weight from the 



FIG. 4: Comparison of the calculated, parameter-free 
LDA+DMFT(QMC) spectra of SrV0 3 (solid line) and 
CaV03 (dashed line) with bulk-sensitive high-resolution PES 
(SrVOa: circles; CaVOa: rectangles) |f| (left panel) and 
Is XAS spectra (SrVO;*: diamonds; Cao.9Sro.1VO3: trian- 
gles) @ (right panel). Horizontal line: experimental subtrac- 
tion of the background intensity. 

quasiparticle peak to the Hubbard bands, and minor dif- 
ferences in the positions of the Hubbard bands. Clearly, 
the two systems are not on the verge of a Mott-Hubbard 
metal-insulator transition. The many-particle DOS of 
the two systems shown in Fig. [21 are seen to be quite sim- 
ilar but not identical. In fact, SrV03 is slightly less cor- 
related than CaV03, in accord with their different LDA 
bandwidth. The inset of Fig. [3] shows that the effect of 
temperature on the spectrum is small for T < 700 K. 
Due to the high photon energy the PES transition ma- 
trix elements will not significantly affect the distribution 
of spectral weight in the energy range considered here. 

In the left panel of Fig.g|we compare our LDA+DMFT 
spectra (300K), which were multiplied with the Fermi 
function at the experimental temperature (20 K) and 
Gauss broadened with the experimental resolution of 
0.1 eV 0, to the experimental PES data obtained by 
subtracting estimated surface and oxygen contributions. 
The quasiparticle peaks in theory and experiment are 
seen to be in very good agreement. In particular, their 
height and width are almost identical for both SrV03 
and CaV03. The difference in the positions of the lower 
Hubbard bands may be partly due to (i) the subtraction 
of the (estimated) oxygen contribution which might also 
remove some 3d spectral weight below —2 cV, and (ii) 
uncertainties in the ab-initio calculation of U. In the 
right panel of Fig. 01 we compare to the XAS data 0] . 
We consider core-hole life time effects by Lorentz broad- 
ening the spectrum with 0.2 eV (23|, multiplying with 
the inverse Fermi function (80K) , and then Gauss broad- 
ening with the experimental resolution of 0.36 eV |24| . 
Again, the overall agreement of the weights and posi- 
tions of the quasiparticle and upper t 2g Hubbard band is 
good, including the tendencies when going from SrV03 to 
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CaVC>3 (Cao.9Sro.1VO3 in the experiment ). For CaVC>3 
the weight of the quasiparticle peak is somewhat lower 
than in the experiment. In contrast to the one-band Hub- 
bard model calculation, our material specific results re- 
produce the strong asymmetry around the Fermi energy 
w.r.t. weights and bandwidths. Our results also give a 
different interpretation of the XAS than in [f| where the 
maximum at about 2.5 eV was attributed to a e g band 
and not to the t 2g upper Hubbard band 0. The slight 
differences in the quasiparticle peaks (see Fig. OJ lead 
to different effective masses, namely m* /m® = 2.1 for 
SrV0 3 and m*/m = 2.4 for CaV0 3 . These theoretical 
values agree with m*/mo = 2 — 3 for SrV03 and CaV03 
as obtained from de Haas-van Alphen experiments and 
thermodynamics 0, E3 ■ We note that the effective mass 
of CaV03 obtained from optical experiments is somewhat 
larger, i.e., m*/ m o = 3.9 26] . 

In summary, we investigated the spectral properties 
of the correlated 3d 1 systems SrV03 and CaV03 within 
the LDA+DMFT(QMC) approach. Constrained LDA 
was used to determine the average Coulomb interaction 
as c7 = 3.55 eV and the exchange coupling as J =1.0 eV. 
With this input we calculated the spectra of the two sys- 
tems in a parameter-free way. Both systems are found 
to be strongly correlated metals, with a transfer of most 
of the spectral weight from the quasiparticle peak to the 
incoherent upper and lower Hubbard bands. Although 
the calculated DMFT spectra of SrV0 3 and CaV0 3 are 
slightly different above the Fermi energy the spectra be- 
low the Fermi energy are the resulting PES are very 
similiar, the quasiparticle parts being almost identical. 
Our calculated spectra agree very well with recent bulk- 
sensitive high-resolution PES || and XAS [|J, i.e., with 
the experimental spectrum below and above the Fermi 
energy. Both compounds are similarly strongly corre- 
lated metals; CaV03 is not on the verge of a Mott- 
Hubbard transition. 

Our results are in striking contrast to previous theo- 
ries and the widespread expectation that the strong lat- 
tice distortion leads to a strong narrowing of the CaV03 
bandwidth and, hence, much stronger correlation effects 
in CaV03. While the c g bands indeed narrow consid- 
erably, the competition between decreasing d—p—d and 
increasing d—d hybridization leads to a rather insignifi- 
cant narrowing of the t2 S bands at the Fermi energy. This 
explains why CaV03 and SrV03 are so similar. With our 
theoretical results confirming the new PES and XAS ex- 
periments, we conclude that the insulating-like behavior 
observed in earlier PES and BIS experiments on CaV03 
must be attributed to surface effects |2flj . 
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